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ABSTRACT 


^/^e^odvnamic developments for vertical axis and horizontal axis wind turbines are given that relate to the 
perfomance and aerodynamic loading of these machines. Included are: (1) a fixed wake aerodynamic model 

of the Darrieus vertical axis wind turbine; (2) experimental results that suggest the existence of a lam- 
inar flow Darrieus vertical axis turbine; (3) a simple aerodynamic model for the turbulent windmil 1/vortex 
ring state of horizontal axis rotors; and (4) a yawing moment of a rigid hub horizontal axis wind turbine 
that is related to blade coning. 


I . DARRIEUS FIXED WAKE THEORY 

A fixed wake theory for the Darrieus Rotor has been 
developed and compared to test results (1). A 
comparison of theory and test results 1s shown 
below in Figure 1. 



The Induced velocities at the front and rear of 
the streamtube shown in Figure 2 are determined 
from the shed vorticity in the wake. For the 
front position (segment BC), a semi-infinite 
wake of strength directly proportional to the 
circulation at BC is developed. No effect from 
the rear position is felt along BC. 

At the rear position, DA, the induced velocity is 
the result of the wake from the front blade which 
appears as an infinite wake and a semi-infinite 
wake of strength directly proportional to the 
circulation at DA. A point in the far wake is 
subject to infinite wakes from both the front 
and fear blade positions. If we denote the in- 
duced velocity by AUp = ap and AUp = 
where F stands for front and R for rear and we 
use Tp and to represent the front and rear 
blade’ circulations, then we may write 



Figure 2. Fixed Wake Streamtube. 


AUp = klFpl = apV„ 

= k[2|rpl + lr^|] = 

= ^[2|rp| + 2lr,|3 = 2av^ 

Here the far-wake induced velocity is defined as 
2a. Accounting for the sign of the circulation, 
these results may be put in the form 
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These relations are exactly the same as developed 
by Holme (2). Momentum considerations yield an 
expression for the dimensionless induced velocity 


an.a) = M rhZ 
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where B is the number of blades, c is the chord, 

X is the tip speed ratio, C. 1$ the lift coeffi- 
cient and W is the velocity relative to the blade. 
Both the lift coefficient and the relative veloci- 
ty are determined fro m t he local induced veloci- 
ties. When the circulation is expressed in terms 
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of the 11ft coefficient, the above three equations 
can be used to determine a, and ap. 

For idealized aerodynamics, = 2Trmsina, the 
above relations can be solved in simple form. 

The resul ts are: 

a(l-a) = ^ COSY sine(l-a + /T^) 

= ^(1 - /I -2a) 

= i(1+2a - /T7a) 

where y is the inclination of a blade element with 
respect to the vertical and a is the blade angular 
position, 90° being directly upwind. 

IK DARRIEUS DRAG COEFFICIENT MEASUREMENTS 

During testing of the 5 meter and 17 meter Darrieus 
vertical axis turbines at Sandia Laboratories , data 
has been obtained at v/1nd speeds below 1 mph (3). 

At such wind speeds, the tip speed ratio is 70 or 
above so that the blades are operating with very 
small angles of attack. In this mode of operation, 
power input is required to overcome blade drag, so 
that the power input, RPM and rotor geometry can 
be used to determine the blade drag coefficient 
Cdq, Data from the Sandia Laboratories test site 
is^reduced using the method of bins (4). The 
results for NACA 0015 blades are shown below. 

Each point represents over 400 samples using the 
bins method. Also shown are the results of a 
numerical prediction for the NACA 0015 airfoil. 

The computation method is described in Reference 
(5). 



REYNOLDS NUMBER 


Figure 3. versus Reynolds Number for the 
' NACA 0015 Airfoil. 

These few data points suggest that laminar f1pw may 
exist up to a chord Reynolds Number of 1.1*10° for 
the NACA 0015 extruded aluminum Darrieus vertical 
axis wind turbine blades. The test results for 
the 17 m at 38.7 RPM show a higher peak power co- 
efficient (0.46) and lower peak power output 
(K = 0.0080) than achieved for operation at 
^max 


higher RPM. If laminar flow characteristics can 
be incorporated into the design of a vertical 
axis wind turbine, significant improvements can 
be made in the COE and reliability of a Darrieus 
Rotor, 

III. TURBULENT WINDMILL/VORTEX RING STATE 

Strip theory analysis of horizontal axis wind 
turbines determines local induced velocity by 
equating blade normal forces to the momentum 
change in a streamtube. Figure 4 below illus- 
trates the operating states of a rotor expressed 
in terms of the thrust coefficient, Cj, and the 
induced velocity, a. 



Figure 4. Rotor Operating States. 

For operations at induced velocities greater than 
0.5, it may be seen that there is a great differ- 
ence in the value of Cj between the momentum and 
Glauert (6). 

The aerodynamic loads and the performance predicted 
using the momentum and Glauert values of Cj show 
corresponding differences. The Glauert empirical 
curve js represented by an implicit algebraic 
equation, however, approximating the Glauert re- 
lation by a straight line gives a useful and ac- 
curate representation the local thrust coefficient 
for strip theory use. Using the expressions 
given below, strip theory calculations were made 
for a two-bladed rotor with NACA 0015 constant 
chord, untwisted blades. The local values of Ct 
used were 

Cj = 4a{1-a) a < 0.38 

' I 


C-P = 0.5776 + 0.96a a > 0.38 

‘l 

Test data (7) and calculations are illustrated in 
Figure 5 below. The measured value of Cj)p 
increases with tip speed ratio for 0° and 
2° blade pitch angle cases. The calculated flow 
for these cases are in the turbulent windmill 
state. Agreement between theory and test suggests 
the above relations can be used to predict the 
aerodynamic loads of horizontal axis wind turbines 
operating in the turbulent windmill state. 


4 





-igure 5. Measured and Predicted Rotor Thrust 
Coefficients at Various Blade Pitch 
Anql es . 

IV. YAW MOMENT FOR A HORIZONTAL AXIS WIND 
TURBINE 


Downwind rigid-hub horizontal axis wind turbines 
dertve their yaw stability from blade coning. 

The coning acts similar to dihedral in an airplane 
by increasing the velocity normal to the upwind 
blade. Yaw forces developed by rigid hub rotors 
are quite small so that the yaw moment controls 
yaw stability. The yaw moment can be expressed 
in terms of Integrals over the blade length. 
Denoting the yaw angle as 6 and the position 
of the 1th blade by 6., as shown in Figure 6, 
the rotor yaw moment *is the moment about the z 
axis. Considering the flow Induced by rotor yaw, 
the yaw moment is obtained by the approach of 
Ribner (8) in which the yaw moment is determined 
from blade forces and momentum considerations. 


WIND 



Figure 6. Wind Turbine in Yaw. 


The results for a rotor with B blades is 

n ? hho 

'■ - ^ - - Z cos 6. 

IT Y 1 1-1+14 


Mv " 1 o 3 
Yaw 


where 


aC-r 


dn 


■- =/„ " 


and = {|)[2C^n-a) + nX] 


and = C|^cos$ + CpSini^ 

Cn = C|^ coS(|) + Cjj sin4> - 
^ a a 


Ct = Cj^s1n(t) - C^cos(^ 

Here the local blade coning angle is ip, the blade 
chord is c and (j> is the angle between the local 
relative wind and the plane of rotation. The 
variable n is the radial position r divided by the 
rotor radius R. For a rotor with 3 or more blades, 
the yaw moment becomes 


B6 




- 2 tt I,+I^ 


The above relation has not been verified by com- 
parison to test data, however, downwind free-yaw 
rigid rotors are observed to be stable in yaw and 
the above term has been found to be, by far, the 
largest stabilizing moment. 

CONCLUSIONS 


The material presented above should be of use in 
the design and design analysis of wind turbines. 
Three of the four developments have been ccsnpared 
with test data. 
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QUESTIONS AND ANSWERS 


R.E, Wilson 


From: T, Cur r in 

Q; Effects of articulation on yaw forces and motions? Order of magnitude with coning? 

A: Professor Miller^s paper^ the Weathervaning of Wind Turbines," (Journal of Energy 

Vol. 3, No. 5, Sept~-Oot 2979, pp , 319-320) , develops expressions for teetering and 

articulated rotors. His paper was the first to note the dihedral like effect of 
coning . 

From: J* Glasgow 

Q: How did you get measurements for the power coefficient on Sandia machine? 

A: Data was obtained using the method of Bins (4) by Sandia (3). 

Fronu F. Perkins 

Q: When will a wind turbine enter. the vortex ring state, and what will the loads look 

like (steady state or not, large or small)? 

A: The blade normal loads as measured (?) are large. The time histories of these 

loads appeared to he steady. 

From: I, Paraschiyoiu 

Qi Did you calculate the loads distribution as a function of the azimuthal angle 9 for 
a vertical-axis turbine? 

Ai Te^, loads were calculated as a function of 6. As would he expected, the upwind 
blade positions experience the largest loads. 

From: L. Mi randy 

Qi Vortex ring formula implies can reach 1.5. X have never seen any data mech. 
over 1.0, have you? 

A: The data of Glauert (British Aeronautical Research Committee Reports & Memoranda 

No. 1026, Feb. 1926) has up to 1.5 and above. 
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